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Abstract-The purpose of the paper is to analyze the damping of 
three-dimensional free vibrations in a transversely isotropic, 
thermoelastic hollow cylinder, which is initially undeformed 
and kept at uniform temperature. The surfaces of the cylinder 
are subjected to stress free and thermally insulated boundary 
conditions. The displacement potential functions have been 
introduced for decoupling the purely shear and longitudinal 
motions in the equations of motion and heat equation. The 
purely transverse wave gets decoupled from rest of the motion 
and is not affected by thermal field. By using the method of 
separation of variables, the system of governing partial 
differential equations is reduced to four second order coupled 
ordinary differential equation in radial coordinate. The matrix 
Frobenius method of extended power series is employed to 
obtain the solution of coupled ordinary differential equations 
along the radial coordinate. In order to illustrate the analytic 
results, the numerical solution of various relations and 
equations are carried out to compute lowest frequency and 
thermoelastic damping factor with MATLAB software 
programming for zinc material. The computer simulated 
results have been presented graphically. 
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I. Introduction 

The vibrations in thermoelastic materials have many 
applications in various fields of science and technology, 
namely aerospace, atomic physics, thermal power plants, 
chemical pipes, pressure vessels, offshore, submarine 
structure, civil engineering structure etc. The hollow cylinders 
are frequently used as structural components and their 
vibrations are obviously important for practical design. The 
investigations of wave propagation in different cylinderical 
structures have been carried out by many researchers [1-6]. 
Ponnusamy [7] studied wave propagation in a generalized 
thermoelastic solid cylinder of arbitrary cross section. Suhubi 
and Erbey [8] investigated longitudinal wave propagation in 
thermoelastic cylinder. Sharma [9] investigated the vibrations 
in a thermoelastic cylindrical panel with voids. 

The theory of elastic vibrations and waves well 
established; see Graff [10] andLove [11]. The objective of the 
present paper is to study the three dimensional vibration 
analysis of simply supported, homogeneous transversely 
isotropic, hollow cylinder of length 'L' and radius 'R'. Three 
displacement potential functions are employed for solving 
the equation of motion and heat equation. The purely 
transverse wavesget decoupled from the rest of motion and 
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are not affected by thermal field. By using the method of 
separation of variables the model of instant vibration problem 
is reduced to four second order coupled ordinary differential 
equations in radial coordinates. One of the standard 
techniques to solve ordinary differential equations with 
variable coefficients is the Frobenius method available in 
literature, see Tomentschger [12]. The secular equation which 
governs the three dimensional vibration of hollow cylinder 
has been derived by using Matrix Frobenius method. The 
numerical solution of secular equation has been carried out 
by MATLAB programming to compute lowest frequency and 
thermoelastic damping factor which have been presented 

graphically with respect to the parameter t, for first two 
modes of vibrations (n = 1,2) . 

H. FORMULATION OF PROBLEM 

We consider a homogeneous transversely isotropic, 
thermal conducting elastic hollow cylinder of length L and 

radius R at uniform temperature T Q in the undisturbed state 

initially. The basic governing equation of motion and heat 
conduction for three-dimensional linear coupled 
homogeneous and transversely isotropic thermoelastic 

cylinder in cylindrical co-ordinates (r , , z) system, in the 
absence of body force and heat source, are given by 
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Here u= (u r , u g , u .) is the displacement vector; 
T(r, 9, z, t) is the temperature change; C u , C l2 , C 13 

and Caa are elastic constants; a,, a 3 and K, , K 3 are 

the coefficients of linear thermal expansion and thermal 
conductivities along and perpendicular to the axis of symmetry 

respectively; p and C e are the mass density and specific 

heat at constant strain respectively, e. the strain tensor ; 

CT.. the stress tensor. The comma notation is used for spatial- 
derivatives and superimposed dot denotes time derivatives. 
We introduce potential functions y/ , G , W as used by 
Sharma [5] 

r r 

Using equation (4) inequation (1) we find that G ,W ,y/ ,T 



satisfies the equations: 
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We consider the free vibrations of a right circular hollow 
cylinder subject to traction free and thermally insulated or 
isothermal boundary conditions on the surface r = R which 

is simply supported on the edges z = an d z — L ■ 
Therefore, we assume solution for three 
displacementfunctions and temperature change as 
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On using solutions (8) in equations (5) and (6), we get 
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Here v 2 is the velocity of purely elastic wave in hollow 

cylinder and <y* is the characteristic frequency. The equation 

(9.1) represents purely transverse wave, which is not affected 
by the temperature change. The equation (9.1) is a Bessel's 
equation and its possible solutions are 
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where jt,' = —A;, • Here E 7 and £"y are two arbitrary 
constants, and / and Y n are Bessel functions for first and 
second kind and / and K n are modified Bessel functions 
for the first and second kind respectively. Generally k , ^ 
, so we go on with our derivation by taking the form of \j/ 

for k, < , the derivation for k, > is obviously similar. 

Therefore the solution valid in case of hollow cylinder is 
taken here as 
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IH. METHOD OFFROBENIUS 

In order to solve the coupled system of differential 
equations (9.2)-(9.4) we shall use matrix Frobenius method. 
The differential equations (9.2), (9.3) and (9.4), can be written 
in matrix form as 

(13) 
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A standard technique tor solving ordinary ditterential 
equations is the method of Frobenius, in which the solutions 

are in the form of power series. Clearly y — is a regular 
singular point of the matrix differential equations (13) and 
hence, we take the solution of type 
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Here the coefficients L k (s),M k (s),N k (s) and the 
parameter s (real or complex) are to be determined. By using 
matrix frobenius method the solution (15) becomes 
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Case-I (when the roots of indicial equation are distinct and 

do not differ by integer) 

Thus the general solution of equation (10) has the form 
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Here E x , E 2 , E 3 , E l , E 2 , E 3 are arbitrary constants 

to be evaluated by using boundary condition. 
Hence 
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Here, potential functions Z(r) written from (8) by using (17a) 

as 
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Here, potential function G, W, X written from (8) by using 
(18a) as 
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Case-II (when the roots of indicial equation are distinct and 

differ by integer) 

Thus the general solution of equation (10) has the form 
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Here E l , E 2 , E 3 , E 4 , E 5 , E b are arbitrary constants 

to be evaluated by using boundary condition.. 
Hence 
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Here, potential functions Z(r) written from (8) by using 
(17b) as 
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Here, potential function G, W, f written from (8) by using 
(31b) as 
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IV. BOUNDARY CONDITIONS 



We consider the free vibration of hollow-cylinder which 
is subjected to two types of boundary condition at lower and 
upper surface (r = R 3 , R 2 ) 
(a) Mechanical condition 
The surface are assumed to be traction free, so that 

*.=*«=**=* ( 21 ) 

(b)Thermal condition 

The surface are assumed to be thermally insulated which 

leads to T = (22) 

V FREQUENCYEQUATION 

In this section we derive secular equation for thermoelastic 
hollow cylinder, subjected to traction free, and thermally 
insulated / isothermal boundary conditions at lower and 
upper surface (r = R 1 , R 2 ) The displacement, temperature 
change and stresses are obtained as 
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prime has been suppressed for convenience. 

Case: /-Using the equation (18a) in (24), we obtain the 

temperature gradient and stresses as 
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Using the boundary conditions (21)-(22) and (25a), we 
get systems of eight simultaneously equations in 

E i ,E i ,E 2 ,E 2 ,E 3 ,E 3 ,E 1 ,E 1 , which willhave non- 
trivial solution if the determinant of their co efficient 
vanishes. This requirement of non trivial solution leads to 
secular equation for hollow cylinder. The secular equation 
are obtained as 
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hollow cylinder 

Case: //-Using the equation (18b) in (24), we obtained the 
temperature gradient and stresses as 
' t. 
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Using the boundary conditions (21)-(22) and (25b), we 
get systems of eight simultaneously equations in 

E,,E 1 ,E,,E A ,E,,E,,E 1 ,E' , which will have non- 
trivial solution if the determinant of their Co efficient 
vanishes. This requirement of non trivial solution leads to 
secular equation for hollow cylinder. The secular equation 
are obtained as 
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VI. NUMERICALRESULTS AND DISCUSSION 

In order to illustrate and verify results obtained in previous 
sections, we present some numerical simulation results. For 
the purpose of numerical computation, we have considered 
zinc-crystal like material whose physical data is given below 
(Dhaliwal and Singh [14]). 

qj = 1.625x10" M« T 2 , q 2 =0.362xlO I] Mb" 1 
q 5 = 0.50S xlO I] xVflf Vjj = 0.627 xlO 11 Mrf ' 

C4 = 0.7T0 xM) u J%i * , p = 175 x 10* Nm' 1 deg" ] 
p i = 5.75 xlO* Mm 1 deg" ] T X, =1.24 xl0 : JJ-jm^ deg" 1 
£, = 1:24 xlO 1 ^" 1 deg" 1 c?* = 5.0lxl0 11 ^ 1 

p = 7.14 xlO-'fcgjjs"-' , T = 296K 

Due to presence of dissipation term in heat conduction 
equation, the frequency equation in general complex 
transcendental equation provides us complex value of 
frequency ( CO ). For fixed value of n and k, the lowest 

frequency (Q ) and dissipation factor (D) are defined as 



n = s> r R 



(£ 



D = a>,R 
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where si, = Re (m) and & . = Im (s> ) - 
The thermoelastic damping factor is given by 
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Re{ &) 
The numerical computation has been carried out for 

n>0 , k>0 with the help ofMATLAB files .The secular 

equation (26) has been expressed in the form of Q = g (fi ) 
and the fixed point iteration numerical technique as outlined 
in Sharma [15] is used to find approximate solution of 
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Q = g(Q) near to the initial guess of the root with tolerance 
(X0 5 )■ We computed lowest frequency and thermoelastic 
damping factor of first two modes of vibrations (n = 1,2) 

for different value of parameter axial wave number t of the 

cylinder. The variations of computer simulated lowest 
frequency and thermoelastic damping factor are with 

parameter t, in respect to first two modes of vibration in 

figures 1 to 2. Fig. 1 represents the variations of lowest 
frequency of first two modes of vibrations, versus axial wave 
number for different values of thickness to mean radius ratio 
of the simply supported hollow cylinder of zinc-crystal like 
material, respectively. It is observed that the lowest frequency 
increases monotonically with axial wave number in both first 
and second modes of vibrations though its magnitude has 
been noticed to be larger in second mode of vibrations as 
compared to that of first mode of vibrations. It is also noticed 
that the magnitude of lowest frequency for each modes of 
vibrationsincreases monotonically with increasing value of 
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Fig 1 . Lowest frequency (Q ) of first and second modes 

(n = 1,2 J versus parameter (t L ) for different values of 
thickness to mean radius ratio. 
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Fig 2. Thermoelastic damping factor [Q ) of first and second 
modes (n = 1,2 J versus parameter (f J for different values of 



thickness to mean radius ratio. 



From fig. 2 it is observed that the thermoelastic damping 
factor of each mode of vibrationsincreases monotonically 
with increasing value of axial wave number for different values 
of thickness to mean radius ratio of the hollow cylinder. The 



trends of variations of thermoelastic damping factor for each 
considered mode of vibrations are almost steady and uniform 
for of thickness to mean radius ratio of the hollow cylinder. It 
is also revealed that the variations of thermoelastic damping 
factor for each mode of vibrations is dispersive in the range 
and remains close to each other for in case. Moreover, there 
exists atleast one value of parameterforeach at which the 
magnitude of thermoelastic damping factor is same for each 
considered mode of vibrations. 

CONCLUSIONS 

The modified Bessel functions and Matrix Frobenius 
method have been successfully used to study the vibrations 
of a homogeneous, transversely isotropic hollow cylinder 
based on three-dimensional thermoelasticity after decoupling 
the equations of motion and heat conduction with the use of 
potential functions. The decoupled purely transverse mode 
is found to be independent of rest of the motion and 
temperature change. The various thermal and mechanical 
parameters have significant effects on the natural frequency, 
thermoelastic damping factor of the hollow cylinder. The 
thermoelastic damping factor increases monotonically with 
axial wave number but decreases with thickness to mean 
radius ratio of the hollow cylinder. 
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